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^ 1 Abstract 

The long-term habitability of Earth-like planets requires low orbital eccentricities. A 

^ secular perturbation from a distant stellar companion is a very important mechanism in 

j> exciting planetary eccentricities, as many of the extrasolar planetary systems are associated 

^ with stellar companions. Although the orbital evolution of an Earth-like planet in a stellar 

\Q binary is well understood, the effect of a binary perturbation to a more realistic system 

O containing additional gas giant planets has been very little studied. Here we provide analytic 

cn 

O 

o 



criteria confirmed by a large ensemble of numerical integrations that identify the initial 

^ orbital parameters leading to eccentric orbits. We show that an extra-solar earth is likely to 

experience a broad range of orbital evolution dictated by the location of a gas-giant planet, 

> necessitating more focused studies on the effect of eccentricity on the potential for life. 

• ^ 

X 

^ 2 Introduction 



2.1 Habitability 

Habitability of a planet is conventionally defined as the capability for liquid water to 
be sustained on a planetary surface. A high (or low) level of stellar flux incident upon a 
planetary surface results in the loss of liquid water through evaporation by the runaway- 
greenhouse effect (or through freezing by global refrigeration) (Kasting et al, 1993). Planets 
with eccentric orbits experience varying levels of stellar flux throughout a year as the distance 
from the star fluctuates during an orbit. The time-averaged stellar flux < F > over the 



1 



orbital period is 



L 
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where L is the luminosity of the host star, and a and e arc the semi-major axis and eccentricity 
of the planet, respectively (Williams and Pollard, 2002). As seen in Eq. 1, the eccentricity 
growth of a planet results in an increase in the time-averaged stellar flux and thus eventually 
leads to unhabitable conditions. Climate models show that the Earth would start to lose 
its surface water with a hypothetical eccentricity of 0.4 around the Sun; if the eccentricity 
exceeds 0.7, the Earth would lose all of its liquid water through the runaway-greenhouse 
effect (Williams and Pollard, 2002). Thus, a perturbation that causes a planet to deviate 
from a circular orbit can greatly impact its long-term climatic stability. This in turn could 
disturb the possible origin, evolution, and prevalence of life on the planet. 

2.2 Kozai Mechanism 

A secular perturbation from a stellar companion is one of the most efficient mechanisms 
in exciting planetary eccentricities to very large values. Unlike our own Solar System, at 
least 20 % of the ~ 240 extrasolar planetary systems detected as of 2007 are members of 
multiple-star systems (Raghavan et ai, 2006; Desidera and Barbicri, 2007; Eggcnbcrgcr 
et ai, 2007). Moreover, the multiplicity among the current sample of planet-hosting stars 
should be higher than 20 % as the photometric searches for stellar companions around known 
planetary systems are still ongoing. 

These stellar companions, despite their large distances from the planetary systems (typi- 
cal binary separations range from ~ 10^ - lO'* AU), can still secularly perturb the planetary 
orbits around the primaries through a unique three-body interaction called the "Kozai mech- 
anism" , in which the cyclic angular momentum exchange between a binary companion and a 
planet results in a large-amplitude eccentricity oscillation of the planet (Kozai, 1962; Holman 
et ai, 1997; Innanen et ai, 1997). The Kozai mechanism takes place when the initial relative 
inclination /pi_bin between the planetary and binary orbits exceeds the critical Kozai angle 
Ikoz = 39.23°. Under such an initial configuration, the planet's eccentricity grows from ~ 
and oscillates with an amplitude that is constant to the lowest order as 



(Holman et al., 1997). For example, an Earth-mass planet at 1 AU would reach the habitabil- 
ity limits of e =0.4 and 0.7 if there is a stellar companion initially inclined by /pi-bin = 45° 
and 56°, respectively. Note that the distribution of /pi- bin in space is most likely isotropic 
because the orbital orientation of binaries with separations greater than ~ 100 AU are not 
expected to be correlated with the invariable plane of the planetary systems around the 
primaries (Hale, 1994; Takeda et al., 2008). 

During Kozai cycles, a stellar companion repeatedly applies a small torque on the plan- 
etary orbit. This torque accumulated over many binary orbits results in the precession of 
the planet's peri center argument and hence the oscillation of the orbital eccentricity. The 
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timescale of the pericenter and eccentricity evolution can be analytically estimated as 
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(Kiseleva et al, 1998; Ford et al, 2000), where the subscripts bin, 0, and 1 refer to the 
binary companion, the primary star, and the planet, respectively. If there is another source 
of perturbation that precesses the planetary orbit on the timescale shorter than tkoz given 
in Eq. |3j then the torque applied from the stellar companion to the planetary orbit secularly 
averages out to zero, and no eccentricity growth by the companion would be observed (Wu 
and Murray, 2003). 

The Kozai timescale may largely vary depending on the system parameters; for example, 
an Earth-size planet with a solar-mass stellar companion at 750 AU would experience Kozai 
oscillations with a period tkoz ~ 320 Myr. If a stellar companion is closer, then the oscillation 
timescale would be largely reduced; for example at 250 AU, tkoz ~ 15 Myr. Note that this 
oscillation period is also relevant to the evolution of life on a planet. For example, if the 
Earth had an eccentricity of 0.4, even though some liquid water would remain on the surface, 
the surface temperatures would exceed 70° C near the pericenter of the orbit (Williams and 
Pollard, 2002). At this temperature, thermophiles and other simple organisms may be able 
to comfortably survive, but complex life could be eliminated (Levy and Miller, 1998; Daniel 
and Cowan, 2000). However, if the eccentricity oscillation period is sufficiently long (on the 
order of hundred million years), the low eccentricity duration of the cycle may provide an 
intriguing avenue for the evolution of more complex organisms. 

Although the eccentricity variation of a planet in a three-body system is understood well, 
in reality, a single Earth-like planet in a binary is probably not a common configuration. 
Both observations and numerical planet-formation simulations rather suggest that planetary 
systems naturally form in multiple configurations, containing one or a few gas giant planets 
(Ida and Lin, 2005; Thommes et al, 2008). Thus, it is certainly relevant to discuss the 
dynamical evolution of a habitable Earth-like planet, secularly perturbed by other planets 
and a distant stellar companion. However, the evolution of gravitationally-coupled multiple 
planets under the infiuence of a secular binary perturbation is significantly more complex 
than the simple three-body Kozai mechanism and has been very little understood in previous 
studies. 



3 Methodology 

3.1 Numerical Simulations 

Here we investigate a prototypical four-body system, containing the Sun, an Earth-mass 
planet, a Jupiter-mass planet, and a hypothetical stellar companion. To determine the 
dependence of the earth's eccentricity on the orbital elements of the jupiter and the stellar 
companion, we have performed an extensive set of numerical simulations of four-body systems 
with different initial binary parameters, /pi-bin and aun, as well as the initial location of the 
jupiter ajup. To ensure that there is no immediate instability in the system, we located the 
jupiter safely within the Hill stable limits for all the ensembles (Gladman, 1993). 
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The orbital evolution of the four-body system was followed by integrating the full equa- 
tions of motion with the Bulirsch-Stoer integration scheme (Chambers, 1999). We have also 
included the post-Newtonian correction to observe the general relativistic (GR) effect on the 
earth that occasionally competes with the secular perturbations from the other bodies. After 
numerically integrating each system for 500 Myr, we have recorded the maximum eccentricity 
attained by the earth during the evolution. 

The post-Newtonian corrections are relevant in these simulations because a wide bi- 
nary companion's perturbation can be weaker than the GR precession of the earth. For a 
Sun/Earth system with a solar-mass binary companion, the GR precession overtakes the 
binary perturbation when the semimajor axis of the companion is greater than ~ 500AU. 
Note that in our simulations this situation does not ensure a stable circular orbit for the 
earth, because of the presence of another planet. 

The 500 Myr simulation is sufficiently long enough to catch all of the systems that would 
result in an eccentric earth. 500 Myr is longer than the characteristic perturbation timescales 
derived from the Kozai mechanism or the Laplace-Lagrange secular theory. For the eccentric 
cases, the maximum eccentricity of the earth was reached quickly within 50 Myr, while the 
stable- circular orbits showed no indication of ever becoming eccentric. 

Each system was integrated 10 times with random initial phases to confirm that the 
initial phase does not affect the overall behavior. The maximum eccentricity of the earth 
that resulted was very similar for each initial phase due to the secular nature of the dynamical 
behaviors. For Figure 1, for each system, the average of the maximum eccentricities of the 
10 simulations was Tised. 

Figure 1 illustrates how the earth's maximum eccentricity depends on ajup and /pi-bin, 
with the binary separation fixed at 750 AU. We first note that without the presence of the 
jupiter, this map would look completely dark; the stellar companion perturbation on the 
earth is too weak - the Kozai oscillations {tkoz ~ 100 Myr) are slow compared to the GR 
precession {tgr ~ 30 Myr) - such that the orbit of the earth would remain nearly circular, 
no matter how large the initial inclination is. If the companion perturbation was stronger 
than GR {tkoz < tgr), without the jupiter, the earth would have experienced eccentricity 
oscillations with amplitudes determined by /pi-bin, above /koz- With the presence of a jupiter, 
however, the earth's orbital evolution is not simply constrained by /pi_bin or the strength 
of the companion perturbation on the earth. Rather, the location of the jupiter and its 
consequent orbital evolution sensitively affects the earth's eccentricity and its potential for 
hfe. 

3.2 Timescale Analysis 

This complex parameter dependence exhibited in Figure 1 can be understood by com- 
paring, on both planets, the strength of the binary perturbation with that of the mutual 
perturbation the planets apply onto each other (hereafter denoted with subscripts "koz" 
and "pp", respectively). Each of these perturbations can be quantified by their precession 
timescales r of the argument of pericenter oj and ascending node fl. The binary perturba- 
tion and planet perturbation timescales were derived from the Kozai perturbation theory 
and the Laplace-Lagrange secular theory, respectively (Takeda et al, 2008). To distin- 
guish between the perturbations on the earth and jupiter, we denote the timescales with 
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subscripts ® and " Jup" , respectively. Figure 2 shows the timescales of the various perturba- 
tions computed over the same simulation space shown in Figure 1. This timescale analysis 
is applicable when the binary orbit is relatively wide compared to the planetary system (to 
apply the Kozai timescale estimate), and when the planets are roughly coplanar (to apply 
the Laplace-Lagrange timescale estimate) (Takeda et al, 2008). 

4 Results 

4.1 Analysis of Simulation Results 

Here we will explore the parameter space of Figure 1 and 2, using the timescale analysis 
to justify the observations in the numerical simulations. 

Figures 3-7 are typical simulation results from each of the dynamical regions to be 
described. Throughout these figures, the lighter marks indicate the Jupiter's parameters, 
while the heavy marks indicates the earth's parameters. The figures show how the planets' 
ascending node fl, argument of pericenter cu, inclination i, and eccentricity e change over 
time, from top to bottom, respectively. The inclination is measured using the binary plane 
as the reference plane. The dotted line in the inclination frame shows the mutual inclination 
between the planets; because the ascending nodes of the planets precess, the difference 
between the inclinations of the planets relative to the binary plane does not indicate the 
planet's mutual relative inclination. 

When the jupiter is sufficiently close to the earth but safely outside the Hill stability 
boundary (1.5 AU < ajup < 3AU; Region A in Figure 1), for both planets, their respective 
planet-planet mutual perturbations dominate over the binary perturbation (rpp jup < Tkoz.Jup 
, Tpp,e < TKoz,e)- ^ result, even if the planets are initially inclined with respect to the 
binary orbit by more than /koz, the Kozai mechanism is completely suppressed, and neither 
planet undergoes eccentricity oscillations in this region. 

When the Jupiter's initial location is moved further from the earth, the Jupiter under- 
goes the Kozai mechanism because the perturbation from the binary companion begins to 
dominate over the perturbation caused by the earth (rKoz,jup < Tpp.jup, Region B and C). 
In Region B, /pi_bin is initially below /koz and thus the orbits of Jupiter and the earth re- 
main nearly circular. In Region C (3 AU < ojup < 13 AU, /pi-bin > 40°), an interesting 
coupling behavior between the earth and jupiter is observed. The Jupiter is subject to the 
Kozai mechanism caused by the stellar companion; its orbital inclination and eccentricity 
oscillate on a timescale Tkoz.Jup- In the simulations, we also observed the earth's eccentricity 
and inclination oscillate in concert with those of the Jupiter, also on the timescale tkoz.Jup- 
This is a result of the gravitational couphng between the planets; specifically, the strength 
of the nodal precession by the planet-planet perturbation (Tn,pp < tq^koz.Jup, Tn,Koz,e) forces 
the two planetary orbits to precess together, while maintaining mutually coplanar orbits. 
Once /pi_bin exceeds 70°, the large orbital eccentricity on the Jupiter induced by the binary 
companion results in the two planetary orbits crossing. Crossing orbits typically scatter the 
earth off its original orbit to an orbit with an eccentricity much beyond the habitable limit. 

When the orbit of the jupiter is further widened to the range 13-25AU, another dy- 
namically interesting Region D arises. Notice that in this region the earth gains significantly 
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large orbital eccentricity, even though initially the binary and the planetary orbits are below 

the critical Kozai angle (/pi-bin < -^Koz)- Because the earth and jupiter are widely separated, 
the planet-planet perturbation is not sufficiently strong to maintain the coplanar precession 
observed in Region B and C (tq,koz,Jup < '^^,pp)- As a result, the earth's nodal precession 
lags behind that of the jupiter, generating a nodal offset angle AQ between the two orbits. 
To compensate, the nodal precession of the earth begins to accelerate, effectively minimizing 
AQ. Nodal acceleration of the earth results in a increase in its orbital inclination to conserve 
the component of the orbital angular momentum normal to the binary orbit. The growth of 
the earth's orbital inclination with respect to the binary orbit also generates a large mutual 
inclination angle between the two planetary orbits, beyond /koz- As a result, the jupiter be- 
gins to apply a Kozai-type perturbation on the earth, thus raising to large values without 
the need of an initially high /pi_bin- Although GR precession can suppress such regions of 
high eccentricity, if the outer planet was a gas giant, this will only occur in extreme cases. 

The nodal precession kicks by the jupiter and the subsequent additional Kozai cycles of 
the earth gradually disappear as the jupiter's perturbation further weakens (ajup > 25 AU). 
As the planet-planet interaction weakens, the earth's orbit becomes stable and circular be- 
cause of GR precession or experiences Kozai cycles induced by the binary perturbation (when 
-^pi-bin > -^Koz)- Rcgiou E in Figure 1 shows an earth with a circular orbit for most of the 
systems. High eccentricities would be seen in this region above /koz, if the GR precession 
did not suppress the binary perturbation. 

4.2 Comparison of Timescale analysis 

Following the timescale analysis, we expect the maximum orbital eccentricities of an earth 
to maintain a similar distribution as the parameters of the outer planet and companion are 
varied. Figure 8 compares the numerical simulations with the timescales for three sets of 
simulations conducted with different binary separations. The numerical boundaries between 
different dynamical classes match well with the analytically estimated boundaries yielded 
from the timescale comparisons. 

4.3 Role of General Relativity 

GR precession competes with all of the perturbations that affect the earth. In Figure 
1, GR has little effect because it is not the dominant interaction throughout most of the 
simulation space. This is due to the small mass of the earth and sun, in addition to the 
relatively large distance between the bodies. Yet, as parameters change, such as the mass of 
the sun, the semimajor axis of the earth, and the mass of the earth, stronger GR precession 
may play a noticeable role in the orbital dynamics of the earth. However, it is difficult to 
consider such situations in the context of habitability because these parameters play a direct 
role in the stellar flux and habitability of the earth, independent of eccentricity. Thus, we 
only considered the possibility of GR suppression in the cases where the planet-planet or 
Kozai interactions become weaker, such as by a smaller, more distant binary companion or 
a gas-giant planet. Figure 9 shows the outer planet parameters at which GR effects can 
suppress the eccentricities observed in the highly eccentric region D on Figure 1. Although 
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the GR precession can theoretically suppress such regions of high eccentricity, if the outer 
planet was a gas giant, this can only occur in extreme cases. 

5 Discussion 

5.1 Using the Timescale Analysis 

The timescale analysis provides a method to quantify the strength of the dynamical 
interactions, thus creating the ability to predict the behaviors that should arise from a 
wide range of planetary systems hosted within a wide-binary The consistency between the 
behaviors seen in the numerical simulations with the behaviors predicted by the timescale 
analysis demonstrates the effectiveness of this approach, as shown in Figure 8. 

Using the timescale analysis, we can show, without exploring the parameter space with 
numerical simulations, that a different range of parameters would result in a similar set 
of orbital behaviors. Figure 10 shows the parameters at which an earth-like planet will 
experience the high eccentricities observed in Region D. The parameters were determined 
by obtaining the specific conditions determined by the timescale analysis. 

5.2 Distribution of Eccentricities 

Most notably, our simulations indicate that an earth is likely to exhibit high eccentricities 
across a wide range of parameters. Figure 11 illustrates the probability of the earth reaching 
the habitable eccentricity limits of 0.4 and 0.7 as a function of the Jupiter's semimajor axis. 
The similarity between the dotted and solid lines at the dynamically rigid and nodal libration 
regions (02 < 20) indicate that the Earth is likely to experience extremely high eccentricity 
oscillations rather than a distribution of high and moderate eccentricities like the single- 
planet case. As shown with the timescale analysis and additional simulation sets, a similar 
distribution should be seen across a range of parameters. Thus, an earth-like planet within 
a binary system will most likely exhibit convincingly non-habitable eccentricities. However, 
it is important to note the length of time at which the planet experiences the high orbital 
eccentricities. As described above, the earth-like planet may experience long periods of 
low eccentricities between its periods of high eccentricities. Such secular behavior must be 
carefully studied in the context of the evolution of life before ruling out the possibility of a 
habitable earth. 

5.3 Conclusion 

Our results demonstrate the variety of orbits that an earth can exhibit within a binary 
system when accompanied by a second planet. With such a high prevalence of binaries, 
in addition to the strong likelihood for multiple planets within a system, we should expect 
eccentric earths even when the companion is distant. As the abihty to search for extraso- 
lar earths increases with observational advances, the possibility of other bodies within the 
system must also be taken into account when discussing habitability. The diverse orbital 
behavior that result from the presence of such bodies highlights the necessity to study in 
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detail extrasolar earths without stable circular orbits. Thus, we call on climatologists, bi- 
ologists, and astrobiologists to carefully consider the broad range of possible eccentricities 
when studying the potential for the origin, evolution, and prevalence of life on a planet. 
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Figure 1: Maximum orbital eccentricity of the earth during 500 Myr of evolution for different 
sets of initial /pi-bin and ajup with fixed initial binary parameters aun = 750 AU and Cbin = 
0.2. For the cases in which the earth collided with the central star, Cmax = 1-0 is recorded. 
The map can be separated into five dynamically-unique regions marked by Region A to E. 
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Figure 2: Secular orbital evolution timescales of the earth and j up iter in a binary as functions 
of ajup. The binary parameters are as in Figure 1. A smaller timescale indicates a faster 
precession, which ultimately implies a stronger perturbation. The Kozai oscillations can 
be suppressed if the planet-planet precession is faster than the Kozai precession. A strong 
planet-planet nodal perturbation allows the planets' orbits to remain coplanar relative to 
each other throughout their evolution. 
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Figure 3: Simulation results of a representative system within Region A, where /pi-bin 
is 70 degrees, and ajup is 3 AU. Although the planets are both inclined over the Kozai 
angle, neither experience eccentricity oscillations because of the strong mutual planet-planet 
interaction. The planets remain coplanar relative to each other. The jupiter is circularized 
by the earth (notice the pericenter precession is ~ Tpp^jup), and the earth is circularized by 
the jupiter. 
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Figure 4: Simulation results of a representative system within Region B, where /pi_bin 

is 20 degrees, and ajup is 9 AU. The planets still remain coplanar, because of the strong 
Tn.pp- The jupiter, inclined below /roz, does not show the eccentricity oscillations although 
it experiences the binary perturbation as seen by its pericenter precession near tkoz.Jup- The 
earth remains circularized by the jupiter. 
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Figure 5: Simulation results depicting the concerted eccentricity described in Region C - 
in this particular system /pi_bm is 50 degrees, and ajup is 9 AU. Notice how throughout 
the eccentricity oscillations, the planets remain coplanar. The jupiter, inclined over /koz, 
experiences the Kozai oscillations. Because of the coplanarity maintained by the strong 
planet-planet interaction, the earth follows the Jupiter's eccentricity oscillations. 



14 




Figure 6: Simulation results of a representative system within Region D, where /pi_bin is 
20 degrees, and ajup is 18 AU. The right plot is the same simulation as the left, but shows 
only the first 20 Myr of the integration. The earth quickly experiences high eccentricities, 
although the planets are both initially inclined below the Kozai angle. As the misalignment 
occurs in the ascending node, the mutual inclination between the planets rises quickly over 
the Kozai angle. This results in high eccentricities induced in the earth, which reappears 
often throughout the 500 Myr. 
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Figure 7: Simulation results of a representative system within Region E, where Jpi^bin is 
50 degrees, and ajup is 31.5 AU. The earth remains at low eccentricity throughout the 500 
Myr because of the GR precession (pericenter precesses ~ 30 Myr, the timescale of the 
GR precession). The jupiter here experiences the Kozai oscillation, because the planets 
are initially inclined over /koz- Although the planets become misaligned, and the mutual 
inclination becomes greater than the Kozai angle, the jupiter is distant enough from the 
earth, such that the GR precession suppress^the planet-planet interaction. 
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Figure 8: The consistency between the timescale analysis and the numerical simulations 
can be seen as the parameters for simulations sets are changed. The three maps on the 
left are numerical simulations of a multiple-planet system with a binary companion, similar 
to the map described in Figure 1. Each map has a different semimajor axis for the binary 
companion. The timescale analyses on the right are derived from the Kozai mechanism and 
the Laplace-Lagrange secular theory - the type of line matches with the type shown in Figure 
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Figure 9: The outer planet parameters at which GR will suppress the Kozai cycles of an 
earth-like planet as described in Region D. The ability for GR to suppress the eccentricities 
observed in Region D will depend on the eccentricity of the outer planet as well. The shaded 
region depicts the largest possible range for the possibihty of suppression. Note how the 
high eccentricities observed in Region D will not be suppressed if the outer planet is at least 
as massive as Jupiter. 
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Figure 10: Calculated semimajor axis and mass of the outer planet for which Region D will 
occur for a system with an earth-like planet. The data is based on calculations derived from 
the timescale analysis for various binary-star separations indicated by the line type. 
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Figure 11: Fraction of systems where the Earth exhibited an eccentricity above 0.4 (dotted) 
and 0.7 (sohd) with respect to ojup, from the simulations performed for Figure 1. 
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